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ABSTRACT 

The  factors  influencing  both  dynamic  hyperinf 
transport  under  the  condition  of  low  volume  high 
pulmonary  ventilation  (HFV) *  were  examined.  In  t 
we  measured  gas  transport  and  dynamic  lung  volumefe  in 
tracheostomized  human  subjects  simultaneously  in  /order  to  help 
determine  the  mechanism (s)  responsible  for  gas  trapping.  We 
found  that  application  of  a  fixed  tidal  volume  (-^-50-100  ml)  over 
a  range  of  frequencies  from  0.5-20  Hz  resulted  in  increased  gas 
transport  efficiency  until  a  critical  frequency  was  reached; 
above  this  frequency  no  further  increase  in  alveolar  ventilation 
was  achieved.  In  contrast,,  mean  lung  volume  (dynamic  functional 
residual  capacity)  increased  only  above  a  critical  frequency. 

These  findings  were  compared  using  a  variety  of  electrical 
analog  models  including  one  which  incorporated  airflow 
limitation;  our  findings  suggest  that  airflow  limitation 
mediates  dynamic  hyper-inflation.  We  also  made  measurements  of 
gas  transport  during  oscillatory  flow  in  a  hardware  model  which 
incorporated  branching  airways. mije  cast  our  results  in  the  form 
of  an  effective  axial  diffusion  coefficient  where 

Deff  =  (K  +  0.5VT1*66f°*94)cm2/sec 

is  the  molecular  diffusivity,  Vm  the  stroke  volume,  and  f  the 
oscillatory  frequency.  At  low  frequency  our  results  are  similar 
to  those  of  Scherer  and  co-workers.  At  higher  frequencies, 
frequency  rather  than  flow  rate  per  se  becomes  an  important 
variable. 


*  abbreviations  and  their  definitions  are  listed  in  the 
glossary  (page  23) 


FOREWORD 


For  the  protection  of  human  subjects  the  investigator (s)  have 
adhered  to  the  policies  of  applicable  Federal  Law  45CFR46. 

Citations  of  commercial  organizations  and  trade  names  in  this 
report:  do  not  constitute  an  official  Department  cf  the  Ar:«y 
endorsement  or  approval  of  the  products  or  services  of  these 
organizations . 


I.  INTRODUCTION 

High  frequency  low  tidal  volume  oscillatory  ventilation 
(HFV)  may  be  of  value  as  a  method  of  respiratory  support  for 
combat  casualties.  However  it  seemed  possible  that 
hyper-inflation  of  the  lung  could  occur  during  application  of 
HFV  and  that  this  hyper-inflation  could  limit  the  utility  of  the 
ventilatory  mode.  Further  since  basic  information  about  the 
mechanisms  controlling  gas  transport  during  HFV  was  not  known, 
studies  were  conducted  in  a  hardware  model  to  better  understand 
how  gas  transport  is  effected  in  such  circumstances. 

Recent  studies  have  demonstrated  that  dynamic  hyperinflation 
may  occur  in  animals  (1)  or  in  patients  (2)  during  low  tidal 
volume  HFV.  In  these  studies,  when  relatively  large  tidal 
volumes  (-100  ml)  were  applied  at  relatively  high  frequencies 
(-15  Hz  or  greater) ,  the  mean  pressure  within  the  alveoli  was 
found  to  exceed  the  mean  airway  opening  pressure,  sometimes  by 
as  much  as  30  cm  H^O.  Simon  et  al  (1)  attributed  this 
phenomenon,  which'  they  observed  in  healthy  dogs,  to  either 
erroneous  pressure  estimation  at  the  airway  opening  due  to 
convective  accelerative  pressure  loss,  or  to  asymmetry  of 
airways  impedance  during  inspiration  and  expiration,  as  a 
potential  explanation  for  this  phenomenon. 

In  cineradiobronchograms  during  HFV  in  dogs  (3) ,  we  have 
noted  nonsinusoidal  diameter  variations  about  the  mean  airway 
diameter  in  a  variety  of  central  airways  during  sinusoidal 
forcing.  In  addition,  pressure  tracings  from  measurement  at  the 
airway  opening  in  patients  (4)  show  similar  asymmetry,  with 
larger  negative  deflections  than  positive  deflections  from  the 
mean.  These  results  strongly  suggest  that  asymmetry  of 
respiratory  impedance  may  be  important  in  creating  disparate 
alveolar  and  airway  opening  mean  pressures. 

The  present  study  in  human  subjects  was  undertaken  to 
address  the  mechanisms  responsible  for  this  asymmetry.  We 
measured  simultaneous  gas  exchange  and  dynamic  hyperinflation  in 
patients,  and  compared  these  results  with  predictions  from  a 
simple  computational  model.  We  found  that  expiratory  flow 
limitation,  rather  than  fixed  asymmetric  or  variable  asymmetric 
resistance,  best  accounts  for  the  experimental  data  reported 
here . 

Since  hyper-inflation  did  not  occur  at  all  frequency  and 
tidal  volume  combinations  studied,  it  seemed  reasonable  to 
further  examine  the  basic  physics  of  gas  transport  when  small 
tidal  volumes  were  used. 

Several  attempts  have  been  made  to  simulate  the  exchange  of 
gases  through  the  respiratory  system  due  to  transport  of  a 
diffusive  nature  (5,6)  and  due  to  direct  alveolar  ventilation 
(7).  Of  these,  diffusive  transport  models  are  most  appropriate 
for  very  small  tidal  volumes.  Although  such  models  are  useful- - 
in  predicting  some  of  the  observed  frequency  and  tidal  volume 
dependence  of  experimental  data,  they  are  deficient  in  that  th :y 
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are  based  merely  on  reasonable  estimates  of  the  actual  rate  of 
dispersion  within  the  lung  rather  than  analyses  or  experiments 
in  which  the  conditions  of  HFV  are  appropriately  modeled. 

Several  groups  have  recently  conducted  experiments  to  help 
elucidate  and  quantify  the  mechanisms  of  transport  under 
conditions  of  small  volume,  high  frequency  oscillation.  Tarbell 
et  al  (8)  and  Azhar  and  Tarbell  (9)  have  measured  the  dispersion 
of  a  bolus  of  liquid  dye  in  a  five-generation  model.  In 
Tarbell' s  experiments,  oscillations  are  imposed  until  the  bolus 
spreads  over  the  entire  network.  The  fluid  is  then  allowed  to 
drain  and  the  distribution  of  dye  is  measured  as  a  funcrion  of 
the  volume  of  fluid  expelled.  Measurements  of  effective  axial 
diffusivity  are  determined  from  the  volume  variance  of  the 
concentration  distribution  and  therefore  represent  an  overall 
network  dispersion  coefficient. 

Some  preliminary  results  have  recently  been  reported  for 
oscillatory  gas  flow  in  bifurcating  systems  of  similar  systems 
of  similar  geometry  (10,11).  While  these  experiments  for  gases 
and  liquids  all  measure  the  rate  at  which  axial  dispersion 
occurs  in  branching  tubes,  there  exist  some  important 
differences  as  to  the  rate  of  molecular  diffusion  in  gas 
transport.  The  results  reported  so  far,  however,  do  not  permit 
an  assessment  of  the  influence  of  molecular  diffusion. 

The  primary  objective  of  these  experiments  was  to  obtain 
correlations  which  can  be  used  in  simulations  of  gas  exchange  in 
the  lung  by  HFV  which  occurs  over  all  airway  branches  (5,7)  .  As 
described  later,  these  results  are  generally  applicable  to 
airways  from  generations  5-13  of  the  Weibel  (12)  morphometric 
model  of  the  human  lung.  In  addition,  however,  these  results, 
in  combination  with  those  of  other  investigators,  help  us  to 
characterize  better  the  flow  and  gas  transport  processes  and  to 
address  questions  relating  to  the  physical  mechanisms  effecting 
gas  transport. 

II .  METHODS 

A.  Experiments  in  Patients 

Four  patients  with  chronic  respiratory  insufficiency  served 
as  subjects  for  the  study;  patients  were  studied  in  the  supine 
position.  Each  required  long-term  mechanical  ventilatory 
support;  their  characteristics  are  provided  in  Table  1. 

Informed  consent  was  obtained  from  each  patient  or  his/her 
guardian . 

Changes  in  lung  volume  were  assessed  with  a  Respitrace 
impedance  plethysmograph.  Respitrace  belts  were  positioned 
around  the  abdomen  and  around  the  chest;  relative  gains  were 
adjusted  so  that  changes  in  the  sura  of  the  two  signals  were 
proportional  to  changes  in  lung  volume.  Calibration  of  the 
Respitrace  volume  signal  was  performed  by  delivering  breaths  of 
brown  volume  from  a  large  syringe. 
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volume 
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The  apparatus  used  to  deliver  HFV  was  similar  to  that 
referred  to  previously  (13) .  The  outflow  rate  of  a  high 
impedance  bias  flow  (mean  rate  30  L/min)  was  servo-controlled  to 
maintain  mean  airway  pressure  at  the  airway  opening  at  a 
predetermined  level  (3.5-20  cm  H20) ,  independent  of  whether  or 
not  oscillations  were  applied.  Airway  opening  pressure  (Pao) 
was  measured  through  a  lateral  side  arm  at  right  angles  to  the 
15  mm  i.d.  HFV  delivery  tube,  using  a  rapidly  responding 
pressure  transducer  (Ailtech  MS-10) ;  mean  Pao  was  obtained  by 
electronically  low-pass  filtering  the  pressure  signal.  Fifty  or 
100  ml  sinusoidal  oscillations  (frequency  0.5-20  Hz)  were 
provided  by  a  sealed  piston  pump  driven  by  a  linear  magnetic 
motor.  Tidal  volume  was  measured  by  integrating  flow  measured 
with  a  number  2  Fleisch  pneumotachograph  coupled  with  a  Validyne 
MP-45  +  50  cm  H20  pressure  transducer,  calibrated  as  described 
previously  (14)  for  each  frequency-tidal  volume  combination  to 
be  employed. 

The  experimental  protocol  was  as  follows:  End-tidal  CCu 
fraction  (Beckman  LB-2)  was  maintained  within  0.5%  (4.8-5.31) 
with  a  conventional  tidal  ventilator,  using  the  patient's  usual 
tidal  volume  at  an  appropriate  frequency.  The  patient's 
tracheostomy  tube  was  then  disconnected  from  the  conventional 
ventilator,  and  the  patient  was  allowed  to  relax  to  functional 
residual  capacity  (FRC)  (Pao=0)  for  2-5  seconds,  until  lung 
volume  (measured  by  the  Respitrace)  was  stable.  The  subject  was 
then  connected  to  the  HFV  apparatus  (with  the  oscillator  turned 
off) ,  and  allowed  to  reach  a  new  stable  lung  volume  determined 
by  the  prechosen  mean  Pao.  Sinusoidal  oscillations  with 
prechosen  frequency  and  tidal  volume  were  then  applied  for  30-45 
seconds.  In  this  interval,  C02  removal  rate  was  measured  as 
bias  flow  rate  times  mean  bias  flow  C02  concentration  (Beckman 
LB-2) ,  and  further  changes  in  lung  volume  were  measured  by  the 
impedance  plethysmograph.  The  subject  was  then  returned  to  the 
conventional  ventilator  (where  FRC  was  shown  to  be  statically 
determined) ,  and  the  end-tidal  C02  again  measured.  The  entire 
sequence  was  repeated  as  HFV  frequency  and  tidal  volume  and  mean 
Pao  were  systematically  varied. 

B.  Computational  Model 

A  simple  electrical  circuit  (Figure  2)  was  used  to  model  the 
presumed  mechanical  behavior  of  the  HFV  pump,  the  bias  flow,  and 
the  central  and  peripheral  portions  of  the  respiratory  tree.  In 
the  analogy,  current  models  volume  flow,  and  voltage  reflects 
pressure,  with  ground  equal  to  atmospheric  pressure.  The  HFV 
pump  is  thus  a  current  source  (Ipump) ;  the  bias  flow  is  modeled 
with  a  voltage  source  (Vmap) ,  whose  value  represents  mean  airway 
pressure,  coupled  to  the  airway  opening  node  by  a  large 
inductance  (Lbf) .  This  configuration  results  in  a  high 
impedance  of  the  "bias  flow"  to  oscillatory  flow. 

The  respiratory  tree  has  two  parallel  branches:  a  central 
airway  "shunt"  compliance  (Caw) ,  which  reflects  potential  radial 
oscillatory  motion  of  the  central  airway  walls,  and  an  axial 


path  from  the  airway  opening  node  to  the  alveolar  node. 
Inspiratory  and  expiratory  limbs  of  the  latter  path  are  isolated 
by  ideal  diodes  to  facilitate  implementation  of  asymmetry;  the 
inspiratory  limb  consists  of  a  fixed  inspiratory  resistance 
(Rinsp)  leading  to  the  alveolar  capacitance  (Calv) .  To  model 
each  of  three  types  of  asymmetric  airways  impedance,  we  varied 
the  structure  of  the  expiratory  limb:  Fixed  asymmetric 
resistance  was  modeled  usina  a  linear  resistor  -  (Rexp) ;  variable 
asymmetric  resistance  was  modeled  with  a  variable  resistance 
whose  value  varied  inversely  with  lung  volume  [Rexp (Palv+30  cm 
H^O) ] ;  and  expiratory  flow  limitation  was  modeled  using  an 
iaealized  PNP  transistor  to  limit  current  at  a  value 
proportional  to  lung  volume  (Palv) ,  in  series  with  a  resistance 
equal  to  Rinsp.  In  each  case,  Palv  equal  to  zero  represents 
FRC ;  Palv  equal  to  -30  cm  H20  models  residual  volume. 

Component  values  were  chosen  to  reflect  relative  physiologic  of 
values . 

Computer  simulations  of  HFV  were  performed  by  an  iterative 
integration  computer  program  (SPICE2,  copyright  Univ.  of  Calif., 
Berkeley).  Tidal  volume  was  set  to  100  ml;  frequencies  of  1-16 
Hz  were  used.  Mean  airway  pressures  of  0,  5,  and  10  cm  H20 
were  studied  for  each  type  of  asymmetry.  Excursion  extremes  of 
Palv  and  Pao  and  mean  Palv  were  estimated  from  this  simulation. 

C.  Gas  Transport  in  Hardware  Models 

The  objective  of  these  tests  was  to  determine  the  dispersion 
coefficient  for  transport  during  conditions  simulating  HFV  in  a 
hardware  model.  For  this  purpose,  we  set  up  a  steady  condition 
in  which  a  trace  gas  is  infused  at  a  constant  rate  at  the 
"alveolar"  end  of  the  network.  Measurements  are  taken  only 
after  a  quasi-steady  state  has  been  reached  throughout  the 
network.  This  state  is  achieved  once  the  cycle-averaged 
concentration  at  each  measurement  location  is  steady,  which,  by 
conservation  of  mass,  indicates  that  the  cycle-averaged  mass 
flow  rate  of  tracer  through  each  generation  is  uniform.  This 
not  only  simulates  the  conditions  of  HFV  under  steady  state 
conditions  (02  and  C02  being  exchanged  at  constant  rates) 
but  also  increases  the  accuracy  of  our  experiments,  since  the 
measurement  can  be  repeated  over  many  oscillation  cycles  and 
averaged  to  eliminate  extraneous  noise. 

The  model  used  in  the  experiments  consists  of  a  system  of 
rigid,  constant  diameter  tubes,  constructed  from  1  cm  nalgene 
"Y"  connectors,  as  shown  in  Figure  5.  Each  bifurcation  has  a 
branching  angle  of  60  degrees  and  each  tube  segment  has  a 
length-to-diameter  ratio  of  3.5:1.  Flow  oscillations  produced 
by  8  (16  in  the  four-generation  experiments)  synchronously 
driven  pistons  and  a  constant  infusion  of  pure  tracer  gas  were 
both  introduced  at  the  "alveolar"  end  of  the  model  so  as  to 
eliminate  any  uncertainties  concerning  the  distribution  of  flow 
to  each  branch.  Studies  (15)  in  a  model  geometrically  similar 
to  curs  have  indicated  that  the  flow  distribution  among  lateral 
and  medial  branches  varies  with  flow  rata  and  can  be  quite 
sensitive  to  small  geometrical  asymmetries. 


The-tracer  concentration  was  measured  simultaneously  at  two 
sites  using  an  infrared  (IR)  absorption  technique.  The 
attenuation  of  monochromatic  radiation  by  a  gas  with 
absorptivity  of  3  is  given  by  Beer’s  law: 


IQ  exp(-$CAx) 


where  I  is  the  attenuated  intensity  of  radiation,  IQ  is  the 
unattenuated  intensity,  C  is  the  average  volume  fraction  of  the 
absorbing  gas  along  a  beam  path  of  length  Ax.  The  experiments 
were  conducted  using  a  helium-neon  laser  operating  at  3.391  urn 
(Jodon  Model  HN-5)  and  methane  as  the  absorbing  tracer  gas 
because  of  the  strong  coincidence  of  the  methane  absorption  band 
(  3  =  9.4  cm-1atm-1)  and  the  laser  wavelength. 

The  laser  beam  was  chopped,  then  split  into  three  beams. 

One  beam  was  directed  toward  a  PbSe  IR  detector  while  the  other 
two  were  passed  through  glass  coverslips  mounted  in  the 
branching  network  at  one  pair  of  the  locations  identified  in 
Figure  5. 

At  measurement  location  (i) ,  the  path  length  averaged 
concentration  can  be  written  as: 


C  -  [k.  -  ln<ei/er)]/2»&  (2) 

where  e^  is  the  output  voltage  of  the  reference  detector,  a  is 
the  light  path  length,  e^  is  the  output  voltage  of  the  optical 
detector  at  measurement  location  (i),  and  (kj_)  is  a  constant 
which  depends  upon  the  alignment  of  the  optical  system. 

We  assume  that  the  path  length  and  cross-sectional  averages 
are  nearly  equal  to  one  another.  This  assumption  breaks  down  in 
situations  in  which  the  cross-sectional  variation  in 
concentration  is  large.  Even  in  these  situations,  however,  we 
might  expect  the  error  at  any  two  of  the  measurement  locations 
identified  in  Figure  5  to  be  in  the  same  direction,  in  which 
case  the  errors  would  tend  to  cancel  when  calculating  the  local 
concentration  gradient. 
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compute  -instantaneous  values  of  the  tracer  concentrations 
(equation  [2]).  The  signal  indicating  piston  position  was  used 
to  compute  the  bulk  oscillatory  flow  rate.  Each  of  these 
signals,  representing  a  time-varying  profile  over  a  1.5  cycle 
period,  was  then  stored  in  digital  form. 

Samples  of  typical  concentration  profiles  at  the  various 
measurement  locations  are  shown  in  Figure  6.  The  subscripts 
refer  to  measurement  locations  identified  in  Figure  5. 

To  conduct  an  experiment,  the  oscillatory  flow  and  tracer 
infusion  were  started  simultaneously,  using  a  methane  flow  rate 
selected  to  achieve  methane  concentrations  of  1-15%  at  the 
measurement  sites.  Concentrations  were  monitored  until  a  steady 
state  was  reached  in  which  the  cycle-average  concentrations  no 
longer  varied  with  time.  At  that  point,  measurements  were 
simultaneously  made  at  two  locations.  When  the  measurement  was 
completed,  the  oscillations  were  stopped,  the  system  was  purged 
with  fresh  air  and  a  reference  measurement  was  taken 
corresponding  to  zero  concentration. 

Computer  analysis  of  the  signals  from  the  IR  detectors 
yields  values  of  methane  concentration  as  a  function  of  time  at 
the  two  measurement  sites.  In  addition,  the  time-mean  methane 
flow  rate  in  each  tube  segment  is  known,  since  a  steady  state 
had  been  reached  prior  to  measurement.  Using  measured  values  of 
the  methane  flow  rate,  q^,  values  for  an  effective  diffusivity 
can  be  computed  from  the" following  expression  (see  reference 
(16)  for  complete  derivation) . 

D.ff  -  ^(Xi-Xjll.Kl-C^/d-Cj))  (3) 


where  subscripts  "1"  and  "2"  refer  to  two  different  measurement 
sites  (either  1  and  2  or  2  and  3  in  Figure  5) .  In  writing  this 
expression,  we  make  the  implicit  assumption  that  the  transport 
process  is  diffusive  in  character,  as  has  been  done  by  others  in 
the  past  (8,17).  However,  neither  our  measurements  nor  those  of 
others  provide,  definitive  evidence  of  the  diffusive  nature  of 
the  transport  process.  In  fact,  Haselton  and  Scherer  (18,19) 
have  demonstrated  that  the  results  of  Scherer  et  al.  (17)  can  be 
explained  on  the  basis  of  a  convective  rather  than  diffusive 
process . 

The  values  of  De£.=  determined  from  Equation  (3)  are  fit, 
by  a  double  least-squares  method,  to  a  correlation  of  the  form: 
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oscillation  frequency.  While  other  forms  might  be  more 
appropriate,  as,  for  example,  that  obtained  by  the  theoretical 
analysis  of  Erdogan  and  Chatwin  (20)  for  steady  flow  dispersion 
in  curved  tubes,  the  correlation  obtained  using  Equation  (4) 
provided  a  satisfactory  fit  of  the  data.  This  form  is 
consistent  with  that  used  by  Tarbell  et  al  (8)  to  reduce  their 
experimental  results  except  that,  for  their  experiments  in 
liquid,  it  was  not  necessary  to  include  the  effects  of  axial 
molecular  diffusion. 

III.  RESULTS 


A.  Experiments  in  Patients 

The  dependence  of  dynamic  hyperinflation  on  frequency  and  on 
mean  Pao  is  illustrated  in  Figure  3.  Using  a  100  ml  sinusoidal 
tidal  volume,  when  mean  Pao  was  low  f(2-5  cm  H20) ,  dynamic 
hyperinflation  occurred  at  the  lowest  frequency  applied  (0.5  Hz) 
and  increased  with  frequency  so  that  a  greater  than  1  liter 
change  in  dynamic  FRC  was  seen  in  two  subjects.  In  contrast, 
when  the  mean  airway  pressure  was  raised  (7-13  cm  H20) ,  a 
different  dependence  of  dynamic  hyperinflation  on  frequency  was 
observed.  Although  the  elevated  mean  Pao  was  itself  associated 
with  an  elevated  mean  lung  volume  compared  to  resting  FRC,  there 
was  no  further  increase  in  lung  volume  as  frequency  was  raised, 
until  the  frequency  reached  some  critical  value  (1.5-4  Hz). 

Above  this  critical  frequency,  lung  volume  did  increase  with 
frequency  in  a  fashion  parallel  to  that  seen  at  the  lower  mean 
Pao . 


In  contrast  to  the  importance  of  mean  alveolar  pressure 
(MAP)  in  determining  dynamic  hyperinflation,  changing  MAP  over  a 
twofold  range  had  little  effect  on  the  dependence  of  rate  of 
CO 2  removal  during  HFV  bursts  on  frequency.  Figure  3  shows 
the  curvilinear  increase  and  plateau  of  VC02  with  frequency 
seen  in  these  patients.  These  results  are  consistent  with  those 
found  in  prior  investigations  (4,13). 

B.  Predictions  of  Transistor  Model 
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The  dependence  of  MAP  on  frequency  and  mean  Pao  predicted 
th  each  of  the  three  asymmetric  impedance  models  is 
lustrated  in  Figures  4a-4c.  In  each  case,  MAP  is  estimated  as 
e  mean  of  the  oscillatory  alveolar  pressure  extremes.  When 
e  expiratory  limb  of  the  respiratory  path  shows  a  fixed 
sistance  whose  magnitude  is  greater  than  that  of  the  Rinsp 
igure  4a) ,  mean  Pao  does  not  affect  the  dependence  of  mean 
Iv  on  frequency  at  fixed  tidal  volume.  In  addition,  for  the 
gher  mean  Pao's,  there  is  no  range  of  frequencies  in  which 
an  Palv  is  indeoendent  of  frequency  at  fixed  tidal  volume, 
edictions  obtained  using  a  variable  Rexp  which  varies 
versely  with  lung  volume  (Figure  4b  are  generally  similar, 
ouch  the 
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variable  Rexp  models  are  qualitatively  different  from  the 
dependence  of  mean  lung  volume  (and  presumably  mean  Palv)  on  HFV 
frequency  and  MAP  observed  in  our  patients.  In  contrast,  the 
dependence  of  mean  Palv  on  frequency  and  Map  predicted  in  the 
expiratory  flow  limitation  model  (Figure  4c)  does  qualitatively 
resemble  that  observed  in  our  patients  (Figure  2) . 

C.  Gas  Transport  Model  Experiments 

The  purpose  of  these  experiments  was  to  study  and  quantify 
gas  transport  during  HFV.  The  range  of  experimental  parameters 
was  therefore  selected  to  overlap,  as  much  as  possible,  the 
range  of  HFV  conditions  found  previously  (4)  to  maintain 
adequate  gas  exchange  in  human  subjects.  In  dimensionless 
terms,  the  experiments  span  a  range  of  Peclet  number  (Pe=2Ua/K 
where  U  is  the  cross-sectional  mean  velocity  amplitude,  a  is  the 
vessel  radius,  and  y  is  the  molecular  diffusivity)  from  5-1100 
and  dimensionless  frequency  S  =a  (  (j/k  'fc  where  to  is  the  angular 
frequency  from  1-12.  In  dimensionless  terms,  this  range  of  Pec 
and  B  corresponds  to  tidal  volumes  ranging  from  2-40  cc  and 
frequencies  from  0.2-13  Hz.  At  a  tidal  volume  of  50  cc  and  a 
frequency  of  15  Hz,  this  corresponds  to  generations  5-13  in  a 
lung  with  the  morphology  given  by  Weibel  (12) .  These  results 
cannot  be  directly  applied  to  the  first  several  generations  both 
because  of  this  constraint  on  Pe  and  also  because  the  turbulence 
generated  by  a  tracheal  catheter  or  glottal  aperture  is  absent 
in  our  tests. 

In  the  calculation  of  De^f  using  Equation  (3) ,  the 
time-averaged  concentration  is  used.  Our  measurement  technique, 
however,  had  sufficient  time  resolution  to  discern  time 
variations  in  concentration  during  the  course  of  an  oscillation 
cycle.  Two  typical  traces  of  concentration  versus  time  are 
shown  in  Figure  6. 

Two  series  of  experiments  were  conducted  that  span  roughly 
the  same  range  of  parameters.  In  the  first  (series  A) , 
measurements  were  made  across  the  first  generation  of  the 
branching  tree,  at  positions  1  and  2  (see  Figure  5) .  In  the 
second  (series  B) ,  concentrations  were  measured  across  the 
second  generation  at  2  and  3.  Since  we  anticipated  differences 
in  the  two  data  sets,  we  performed  the  curve-fitting  analysis  on 
each  set  separately.  In  addition,  we  found  that  the  results 
from  series  A  satisfied  different  relations  depending  on  the 
range  of  dimensionless  frequency.  The  correlations  of  all  our 
data  are  presented  in  Table  2,  including  the  single  expression 
obtained  when  all  the  data  from  both  generations  were  analyzed 
together . 

Shown  also  in  Table  2  are  values  for  a  second  set  of 
constants  corresponding  to  the  dimensionless  relationship  (see 
dimensional  analysis  below) : 


(Deff-T )  tX  =  q  Pcrf3s 


Basis 
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While-it  is  convenient  to  express  our  results  in 
dimensionless  form,  the  dependence  on  molecular  diffusivity 
suggested  by  this  equation  has  not  been  tested  in  our 
experiments,  since  we  used  a  single  gas  mixture.  Consequently, 
the  dimensionless  transport  coefficient  given  above  may  also 
exhibit  a  dependence  on  Schmidt  number  not  evident  in  our 
experiments . 

The  fit  of  our  data  to^this  particular  form  can  be  seen  by 
plotting  (De£f/  < -  l)/(q  s5  )  vs.  Pe.  According  to  Equation 
(5) ,  this  plot:  has  a  slope  of  r  on  a  log-log  plot  as  shown  in 
Figure  7. 


IV.  DISCUSSION 


A.  Patient  Experiments 

The  data  we  report  here  confirm  and  extend  previous 
observations  (1,2).  In  our  patients,  the  rate  of  alveolar 
ventilation  during  HFV  with  either  50  or  100  ml  tidal  volumes 
increased  with  frequency  until  a  critical  frequency  was  reached; 
above  this  frequency,  no  further  increase  in  tfCQ2  (reflecting 
alveolar  ventilation)  was  found  with  increases  m  oscillatory 
frequency.  In  addition,  there  was  little  change  of  ^Cq2  with  • 
mean  Pao  over  the  approximately  twofold  rane  studied.  Each  of 
these  results  coincides  with  previous  observations  on 
HFV-mediated  gas  exchange  (4,14).  Our  findings  of  asymmetric 
Pao  fluctuations  around  the  mean  and  of  impressive  dynamic 
pulmonary  hyperinflation  also  confirm  previous  reports  in  humans 
(2)  and  in  animals  (1) . 

Two  features  of  our  data  represent  new  observations.  First, 
while  alveolar  ventilation  appears  to  depend  on  oscillatory 
frequency,  dynamic  hyperinflation  is  more  directly  related  to 
expiratory  flow.  In  the  two  patients  in  whom  tidal  volumes  of 
50  and  100  ml  were  applied  using  similar  mean  Pao ' s  the  plateau 
of  VCq2  occurred  at  about  the  same  frequency  for  either  tidal 
volume.  However,  the  frequency  at  which  equal  degrees  of 
dynamic  hyperinflation  occurred  was  about  twice  as  large  when 
the  smaller  tidal  volume  was  used  than  when  the  larger  tidal 
volume  was  applied.  This  finding  indicates  that  it  is  flow, 
rather  than  frequency,  which  determines  dynamic  hyperinflation 
at  a  given  mean  Pao.  Although  previous  investigations  (2) 
suggested  this  result,  the  simultaneous  VCq2  and  lung  volume 
measurements  reported  here  strengthen  this  conclusion.  In 
addition,  cur  finding  that  at  the  same  oscillatory  frequency  and 
tidal  volume  and  similar  mean  airway  pressure,  shorter 
expiratory  times  markedly  enhance  lung  hyperinflation  further 
suggests  that  it  is  expiratory  flow,  rather  than  average  or 
inspiratory  flew,  which  is  related  most  closely  to  this 


Second,  in  contrast  to  its  minor  influence  on  alveolar 
ventilation,  mean  Pao  is  an  important  determinant  of  the 
relationship  between  dynamic  hyperinflation  and  oscillatory  flow 
(Figure  3).  When  mean  Pao  was  higher,  there  was,  of  course, 
static  elevation  of  lung  volume  above  the  resting  FRC.  However, 
further  dynamic  hyperinflation  did  not  occur  at  many  flow  rates 
that  did  increase  lung  volume  above  the  lower  static  lung 
volumes  associated  with  lower  mean  Pao's.  This  finding  suggests 
that  if  asymmetry  of  inspiratory  and  expiratory  impedances 
accounts  for  dynamic  hyperinflation,  then  the  asymmetry  first 
appears  with  increasing  flow  at  relatively  lower  flow  rates  when 
the  lung  volume  is  lower,  and  at  rather  higher  flow  rates  when 
the  lung  volume  is  higher.  Maximal  expiratory  flow  limitation 
behaves  in  this  manner.  These  two  new  observations  suggest  that 
maximal  expiratory  flow  limitation  occurs  during  HFV,  and 
substantially  accounts  for  dynamic  pulmonary  hyperinflation. 

To  examine  other  potential  mechanisms  of  asymmetric 
oscillatory  impedance  during  HFV,  we  applied  the  simple 
electrical  analog  of  mechanical  events  during  HFV  illustrated  in 
Figure  2.  We  considered  three  forms  of  asymmetry  which  might  be 
extant:  fixed  asymmetric  resistance,  variable  asymmetric 

resistance,  and  expiratory  flow  limitation  with  equal 
inspiratory  and  expiratory  resistances.  By  examining 
predictions  from  these  models,  it  is  clear  that  only  expiratory 
flow  limitation  results  in  a  dependence  of  dynamic 
hyperinflation  on  oscillatory  flow  and  mean  airway  opening 
pressure  (assuming  that  mean  alveolar  pressure  reflects  lung 
volume  [Figure  4])  which  resembles  that  found  in  our  patients. 

It  is  of  particular  note,  though,  that  the  other  feature  which 
suggests  asymmetry  of  inspiratory  and  expiratory  impedances, 
namely,  asymmetric  fluctuations  about  the  mean  airway  opening 
pressure,  is  qualitatively  similar  among  the  three  models,  at 
least  for  the  higher  flow  rates  at  which  dynamic  hyperinflation 
occurs.  Therefore,  this  finding  of  asymmetric  airway  opening 
pressure  excursions  alone  does  not  allow  us  to  distinguish  among 
the  mechanisms  potentially  responsible  for  dynamic 
hyperinflation. 

Another  feature  similar  among  the  three  models  is  the 
predicted  distribution  of  flow  introduced  into  the  trachea 
(airway  opening  node) .  In  each  case,  for  a  given  tracheal  tidal 
volume,  the  central  airway  wall  motion  "shunt"  volume  increases 
with  frequency  and  varies  only  little  with  mean  airway  opening 
pressure.  Conversely,  the  peripherally  delivered  tidal  volume 
(VTalv)  falls  with  frequency  above  the  critical  frequency;  thus, 
peripherally  delivered  flow  (frequency  times  VTalv)  rises  with 
frequency  until  it  starts  to  plateau  at  the  critical  frequency. 
Like  central  airway  "shunt"  volume,  peripheral  flow  depends 
little  on  mean  airway  opening  pressure.  The  similarity  in  these 
dependences  of  predicted  peripheral  flow  and  of  observed  Vc~,2 
on  HFV  frequency  and  mean  airway  opening  pressure  at  constant 
tracheal  tidal  volume  is  striking,  but  does  not  necessarily 
imply  a  causal  relationship.  Nonetheless,  that  each  of  the 
three  asymmetric  impedance  models  predicts  similar  dependences 
of  peripheral  flow  on  frequency  and  MAP  means  that  this  behaves:- 


does  not- discriminate  among  these  potential  modes  of  asymmetry. 


Our  data  clearly  indicate  that  flow  limitation  is  most 
likely  responsible  for  dynamic  hyperinflation  during  HFV ;  what 
is  still  not  established  is  the  relationship  of  this  limitation 
to  gas  transport  during  HFV.  To  gain  better  understanding  of 
the  latter  problem,  we  studied  gas  transport  in  hardware  models. 

3.  Hardware  Model  Studies 

1.  Effects  of  Model  Geometry 

One  objective  of  these  experiments  was  to  determine 
expressions  that  can  be  used  in  model  simulations  of  gas 
exchange  during  typical  conditions  of  HFV.  The  extent  to  which 
these  results  are  useful  for  this  purpose  depends  upon  the 
realism  incorporated  into  the  model.  In  our  experiments  several 
geometrical  features  of  the  model  differ  from  those  typically 
found  in  the  lung.  For  example,  the  transition  from  one 
generation  to  the  next  is  abrupt  in  the  model  in  that  the 

corners  are  sharp  and  the  change  in  area  (a  factor  of  two  at 

each  branch)  is  larger  than  one  typically  finds  in  the  airways. 
These  differences  are  likely  to  promote  flow  separation, 
especially  during  inspiration.  The  velocity  fluctuations  and 
eddy  motion  associated  with  the  separated  region  could  influence 
the  rata  of  gas  exchange.  The  same  factors  could  cause  an 
earlier  transition  to  turbulence  in  the  model  than  in  the  lung, 

but  since  the  peak  values  of  Stokes  layer  Reynolds  number  (Re/d 

dimensionless  frequency)  in  our  tests  were  always  less  than  200, 
it  seems  unlikely  that  we  entered  the  turbulent  regime  (2) . 
Furthermore,  we  found  no  evidence  in  our  experimental  results  to 
suggest  a  transition  to  turbulence  as  observed  by  Azhar  and 
Tarbell  (9)  . 

2.  Boundary  Effects 

In  experiments  conducted  in  networks  with  relatively  few 
generations,  one  must  also  consider  the  importance  of  end 
effects.  We  would  expect,  for  example,  that  fluid  near  the  ends 
would  exhibit  less  mixing  due  to  secondary  flow,  since  that 
fluid  resides  in  a  straight  tube  segment  during  a  larger  part  of 
the  cycle  than  fluid  in  the  more  central  regions  of  the 
model.  This  effect  would  be  most  evident  in  measurements  made 
at  the  first  generation,  where  the  mean  fluid  displacement  is 
greatest.  We  found,  however,  that  the  differences  between  the 
first  and  second  generations  (series  A  and  B)  were  relatively 
small.  Although  the  curve-fitting  procedure  yields  somewhat 
different  results  for  the  two  test  series  (Table  2),  when  the 
two  data  sens  are  combined,  there  is  a  relatively  small  decrease 
in  the  accuracy  of  the  curve  fit. 

analyzing  the  data,  we  obtained  the  best  fit  to  the  data 
ding  the  results  from  series  A  into  two  "groups,  depend  in  v 
her  .  was  >  or  <  4.  However ,  we  found  no  clear  advantage 
e  series  b  _  welts  were  processed  in  the  same  way .  Whilz 
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this  seems  to  suggest  that  the  transition  to  quasi-steady  flow 
occurs  at  lower  values  of  S  ,  this  may  simply  reflect  the 
relatively  large  experimental  error  in  these  low  frequency 
measurements  as  well  as  the  narrow  range  of  3  tested.  To 
resolve  these  questions  will  require  more  extensive  tests  at  yet 
lower  values  of  3  . 

3.  Dimensional  Analyses  of  Unsteady  Dispersion 

To  provide  a  framework  for  the  consideration  of  the  effects 
of  unsteadiness  and  molecular  diffusion  to  follow,  we  first  cast 
the  axial  dispersion  coefficient  in  non-dimensional  form.  In 
dimensional  terms,  we  express  (Dpff-<  )  (representing 
transport  above  and  beyond  that  aue  to  axial  diffusion  alone)  as 
a  function  of  six  dimensional  parameters: 

Deff  *  fl(U*  *»  T*'1'  (6) 


where  T  is  the  period  of  oscillation  (T  =  2  Vco  )  ,  and  g 
symbolically  represents  all  the  parameters  associated  with 
airway  geometry.  A  dimensionless  form  of  (6)  is: 


(Def  f--*)  /U^T  -  f2(Pe,  P,  Sc.  G)  (7] 

where  Sc  is  the  Schmidt  number  (Sc  =  V//|<  )  and  G  denotes  all 

dimensionless  parameters  having  to  do  exclusively  with  system 
geometry. 

The  dimensionless  time  in  this  expression  is  the  square  of 
the  ratio  formed  by  the  radial  diffusion  time  (  T ~ad  =  a2/^) 
and  the  cycle  period  (T) .  Therefore  we  can  rewrite  (7)  as: 

(Deff-*)/U2T  -  fn(Pe,  Trad/T,  Sc,  G)  (8) 


When  the  theoretical  results  for  laminar  oscillatory  flow 
dispersion  in  a  straight  tube  (22)  are  expressed  in  the  form  of 
Equation  (8),  we  obtain  the  curve  plotted  in  Figure  8.  The 
comparison  in  Figure  8  provides  information  on  the  importance  of 
G  in  Equation  (8).  Note  that  when  '  rad  alone  is  decreased, 
the  rate  of  dispersion  will  either  decrease  or  increase, 
depending  on  whether  Tracj  is  <  or  >  T,  respectively. 

An  alternate  but  entirely  equivalent  dimensionless  form  can 
be  derived  which  is  useful  in  comparing  our  results  to  those  of 
other  investigators. 


=  (Def  f-V) /tf  =  fn(Pe,  p,  Sc,  G) 


(9) 


Under  certain  circumstances,  one  or  more  of  the 
dimensionless  groups  in  Equation  (9)  can  be  eliminated.  Several 
examples  are  of  particular  interest  in  the  following 
discussions:  (1)  If  the  flow  were  quasi-steady,  then  D  would 

be  independent  of  -  .  (2)  If  the  flow  were  highly  turbulent, 

convective  mixing  would  overwhelm  the  effects  of  molecular 
diffusion,  Sc  would  consequently  have  little  effect  and  <  should 
be  replaced  by  '•>  .  (3)  If  the  flow  were  both  quasi-steady  and 

highly  turbulent,  we  would  expect  D  to  be  independent  of  both  3 
and  Sc. 

4.  Comparison  to  Other  Experiments. 

Because  of  substantial  differences  between  our  experiments 
and  those  of  Tarbell  et  al.  (8)  and  Scherer  et  al.  (17) 
(hereafter  referred  to  as  Tarbell  and  Scherer) ,  direct 
comparisons  are  difficult.  The  major  difference  is  that  while 
our  measurements  are  local  in  the  sense  that  we  determine  Deff 
from  the  concentration  gradient  across  a  single  generation,  Both 
Tarbell  and  Scherer  measured  the  rate  of  dispersion  over  the 
entire  branching  network,  and  their  results  are  therefore 
influenced  to  varying  degrees  by  the  local  transport  process  in 
each  individual  generation. 

Despite  these  differences,  some  useful  comparisons  can  be 
made.  In  order  tc  cast  all  the  results  in  comparable  form,  we 
first  note  that  the  dispersion  expressions  obtained  by  Tarbell 
and  Scherer  represent  a  network  average.  We  will  assume  that 
the  dependence  they  observe  on  Reynolds  or  Peclet  number  (and  a 
for  Tarbell)  is  consistent  throughout  all  generations. 
Accordingly,  the  rate  of  dispersion  is  greatest  in  the  lower 
generations,  where  the  highest  Reynolds  numbers  are 
encountered.  The  network  dispersion  coefficient  that  they 
define,  therefore,  can  be  thought  of  as  a  weighted  average  of 
the  local  values  of  Def;f.  When  casting  their  results  in 
dimensionless  form,  Taroell  uses  the  total  cross-sectional  area 
of  the  second  generation  as  the  approximate  "network  area."  In 
the  following  comparison,  we  assume  that  the  network  dispersion 
coefficients  defined  by  Scherer  and  Tarbell  are  roughly 
equivalent  to  the  dispersion  coefficient  in  their  second 
generations.  While  this  is  somewhat  arbitrary,  the  error  in  the 
approximation  is  likely  to  be  relatively  small  and  will  affect 
only  the  numerical  constant  (e.g.,  q  and  q  in  Equation  [5])  and 
not  the  exponents.  Based  on  these  assumptions,  we  will  express 
the  results  of  Scherer  and  Tarbell  in  terms  of  'the  Peclet  number 
in  the  second  generation  (based  on  tube  diameter) .  We  relate 
Tarbell 's  Reynolds  number  (Rem  =  UQa/  where  UQ  is  the 
velocity  amplitude  in  the  parent  tube  of  the  network)  to  a  local 
one,  using: 


Re  =  ReT/2 


Dimensionless  dispersion  coefficients  are  thus  obtained  for  the 
results  of  Tarbell: 


0.0035  Sc (Pe  Sc)1 ‘82a~° *86 


and  Scherer: 


(10) 


oSs  -  0.9  Pe 


(11) 


While  we  have  introduced  Sc  to  Equation  (10)  to  be 
consistent  with  Tarbell ' s  dimensionless  formula,  in  the 
following  comparisons  Sc  is  set  equal  to  unity. 

Within  the  range  of  Reynolds  number  and  dimensionless 
frequency  investigated  in  both  unsteady  experiments  (which  also 
overlap  the  range  of  Peclet  number  tested  by  Scherer) ,  the 
comparison  of  Figure  9  can  be  made  for  a  gas  mixture  with  Sc  = 

1.  Note  that  in  placing  Tarbell 's  correlations  on  this  graph, 
we  have  implicitly  assumed  that  mixing  by  molecular  diffusion  is 
unimportant  and,  therefore,  that  Tarbell 's  results  are 
unaffected  by  Schmidt  number.  The  validity  of  this  assumption 
is  considered  later. 

More  recently,  Azhar  and  Tarbell  (9) ,  using  methods  similar 
to  those  of  Tarbell,  have  shown  that  their  network  dispersion 
coefficient  DN  takes  on  one  of  two  forms: 

Dn  a  Re1,43^1*64  (12) 


for  low  Re,  and: 


Dn  «  Re1,57  (13> 


for  high  Re.  They  found  that  the  transition  value  of  Re 
depended  on  oscillation  volume  and  ranged  from  2,000  to  9,500  in 
their  experiments. 

5.  The  Effect  of  Unsteadiness 

As  noted  earlier,  the  transport  process  is  quasi-steady 
when  -  racj/T  >  1,  leading  to  a  result  in  which  D^ff  is 
found  to  be  independent  of  .  Our  expressions  tor  Deff  show 
a  clearly  significant  dependence  upon  (or  :  ) ,  at  least  for 

values  of  >4  (see  Table  2) ,  suggesting  that  generally 

-,.,^/T  <  1.  The  data  for  <  4  are  inconclusive  in  this  regard-.. 
largely  because  of  a  limited  amount  of  data  in  this  range.  By 


comparison,  the  results  of  Tarbell  given  in  Equation  (11)  yield 
a  smaller,  yet  non-zero  exponent  for  8  . 

However,  the  more  recent  data  from  Azhar  and  Tarbell  show 
two  distinct  regimes  depending  on  whether  Re  is  above  or  below 
some  transition  value.  Their  high  Re  correlation  is  independent 
of  2  and  therefore  quasi-steady.  The  reason  for  this, 
presumably,  is  that  turbulent  mixing  decreases  the  radial 
equilibration  time  so  that  T  racj/T  becomes  <  1.  The  low  Re 
correlation  from  Azhar  and  Tarbell  bears  a  striking  resemblance 
to  ours  (Table  2)  in  terms  of  the  exponents  on  Re  and  ^  .  In 
both  cases  the  result  appears  to  be  intrinsically  unsteady. 

We  can  also  compare  our  low  frequency  results  to  those  of 
Scherer  et  al .  (1975)  for  steady  flow  through  a  branching 

network.  Taking  all  our  data  from  both  generations  for  3  <4, 
we  find  that  the  ratio  D/  Dg  has  a  mean  value  of  0.50  with  a 
standard  deviation  of  0.30  for  the  24  experimental  points  within 
this  range  of  8  .  It  has  been  suggested  by  Pedley  (22)  that  a 
more  appropriate  unsteadiness  parameter  in  a  branching  network 
is  e  =  uiL/U  (where  L  is  the  generation  length)  .  Choosing  those 
data  for  which  £<  I  (37  data  points),  the  mean  value  for  D/  Ds 
is  0.62  ±  0.22.  Thus  both  the  approximate  agreement  in  form 
between  Scherer's  findings  and  our  low  frequency  results  and  the 
good  numerical  agreement  indicated  by  the  values  for  D/  Ds 
given  above  support  our  claim  that  these  low  frequency 
experiments  are  quasi-steady. 

6.  The  Effect  of  Molecular  Diffusion 

Based  on  the  comparison  in  Figure  9,  there  appears  to  be 
reasonably  close  agreement  between  our  results  and  those  of 
Tarbell  in  view  of  the  approximations  on  which  the  comparison  is 
based.  Even  closer  correspondence  is  found  between  our  data  and 
the  low  Re  data  of  Azhar  and  Tarbell.  All  this  suggests  that 
molecular  diffusion  plays  at  most  a  small  role  in  the  transport 
process  of  these  experiments.  This  conclusion  is  consistent 
with  the  results  of  Knopp  et  al.  (23),  who  found  no  significant 
effect  of  molecular  diffusivity  in  the  rate  of  tracer  clearance 
from  dog  lungs  during  high  frequency  oscillation. 

There  is  reason  to  believe,  however,  that  seme  differences 
might  exist  between  liquids  and  gases,  based  on  experimental 
results  of  dispersion  under  steady  flow  conditions,  even  when 
the  flow  is  turbulent  or  in  the  presence  of  secondary  flow.  In 
a  comparison  between  the  theoretical  prediction  of  Taylor  (24) 
fer  axial  dispersion  under  conditions  of  steady  flow  in  a 
straight  tube  and  experimental  results,  Levenspiel  (25)  has 
attributed  the  observed  differences  to  the  effects  of  molecular 
diffusion.  Gas  and  liquid  experiments  conducted  in  curved  tubes 
with  secondary  flows  also  exhibit  large  differences  in  axial 
dispersion  in  the  range  Re  <  10"*  (van  Anael  et  al,  26)  .  Thus, 
while  the  results  obtained  in  branching  tubes  are  not,  by 
themselves,  sufficient  to  shew  a  significant  effect  cf  Sc, 
experiments  in  other  configurations  suggest  that  some 
differences  may  exist  between  liquids  and  gases. 


7.  Conclusion 


Correlations  have  been  present  that  are  applicable  to  models 
of  HFV  in  the  lung.  These  demonstrate  that  tidal  volume  exerts 
a  greater  influence  than  frequency  of  oscillation,  consistent 
with  the  measurements  of  gas  exchange  during  HFV  in  test 
animals . 

Our  results  show  a  significant  dependence  on  dimensionless 
frequency,  demonstrating  that,  at  least  for  Reynolds  numbers  up 
to  about  1000  and  values  of  3  >4,  the  effects  of  unsteadiness 

are  important.  Our  data  lie  remarkably  close  to  the  results  for 
oscillatory  flow  in  a  straight  tube  but  are  higher.  This 
enhancement  of  axial  transport  can  be  attributed  to  one  or  more 
of  the  following  mechanisms:  increased  cross-stream  mixing,  a 
less  uniform  axial  velocity  profile  (both  of  which  enhance 
augmented  dispersion) ,  and  convective  streaming. 

There  remain,  however,  a  number  of  uncertainties.  Among 
these  are  questions  pertaining  to  the  importance  of  molecular 
diffusion  at  both  low  and  high  Reynolds  numbers  and  the  effect 
of  and  conditions  for  a  transition  from  laminar  to  turbulent 
flow.  In  addition,  the  relative  importance  of  streaming  versus 
augmented  dispersion  in  all  the  experiments  conducted  in 
branching  networks  needs  to  be  examined  further. 


V,  MILITARY  SIGNIFICANCE 

Ventilatory  assist  by  HFV  may  be  used  in  battlefield 
resuscitation  of  chemical  casualties.  Our  studies  have  provided 
basic  information  about  gas  trapping  and  gas  transport  in  this 
circumstance.  This  could  be  of  great  value  in  the* design  of 
practical  portable  ventilators. 
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English 

a  =  Distance  along  a  path 

c  =  Volume  fraction  of  absorbing  gas 

Calv  =  Alveolar  compliance 

Caw  =  Central  airway  shunt  compliance 

cm  =  Centimeter 

Deff  =  Effective  diffusivity 

Dr.  =  Network  dispersion  coefficient 

Ds  =  Diffusion  coefficients  given  by  Scherer  (17) 

Dt  =  Diffusion  coefficients  given  by  Tarbell  (8) 

ei  =  Output  voltage  of  detector  -  experimental 

er  =  Output  voltage  of  detector  -  reference 

f  =  Frequency 

FRC  =  Functional  residual  capacity 
G  =  Geometric  coefficients 
HFV  =  High  frequency  ventilation 
Hz  =  Hertz  (cycles  per  second) 

I  =  Laser  intensity  -  experimental 
Ibf  =  Bias  flow  inductance 
Io  =  Laser  intensity  -  reference 
Ipump  =  Current  source  for  HFV 

Ki  =  Optical  system  constant 
In  =  Natural  logarithm 
MAP  =  Mean  airway  pressure 

Falv  =  Alveolar  pressure 

Pao  =  Airway  opening  pressure 
Pe  =  Peclet  Number 
Re  =  Reynolds  Number 
Rexp  =  Expiratory  resistance 

Rir.sp  =  Inspiratory  resistance 

5a  =  Schmidt  Number 
I  =  Oscillaricn  period 
=  Velocity  (mean) 

V 122  =  Rate  of  CO-  removal  by  HFV 

Vmap  =  Mean  airway  opening  voltage 

V  ~  T  i  q  fi  I  v  c  1  u  m  o 


Mil 
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Glossary  (Continued) 


Greek 


Molecular  diffusivity 
Gas  absorbtivity 
Regression  constant 
Regression  constant 
Regression  constant 
Regression  constant 
Regression  constant 
Regression  constant 

Dimensionless  frequency  based  on  molecular  diffusivity 
Dimensionless  frequency  based  on  gas  viscosity 
Angular  frequency 
Radial  diffusion  time 

Unsteadiness  parameter  (  u  x  length/V) 

Kinematic  Viscosity 


Figure  1 


C02  removal  rate  (VC02  versus  oscillatory  frequency 

during  short  bursts  of  high  frequency  oscillation  at 

fixed  tidal  volume,  from  each  of  three  patienus. 

Triangles  connected  by  solid  lines  present  data 

gathered  with  low  mean  airway  pressure  (-  3  cm  H20) ; 

squares  connected  by  broken  lines  present  data  gathered 

with  high  mean  airway  pressure  (~  10  cm  H20) .  Vc02 

increases  with  oscillatory  frequency  up  to  a  critical 

frequency  (1-2  Hz),  above  which  there  is  no  further 
« 

increase  in  VC02.  There  is  no  consistent  dependence 
of  the  pattern  of  Vco2  versus  frequency  on  mean 
airway  pressure. 
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Figure  2 

Three  models  of  the  mechanical  events  potentially 
responsible  for  dynamic  hyperinflation  during  HFV.  In 
each,  the  mean  airway  opening  pressure  is  maintained  by 
a. pressure  source  (Vmap)  connected  to  the  airway 
opening  node  by  an  inertance,  which  allows  steady 
pressure  to  be  transmitted  while  excluding  oscillations 
from  entering  Vmap.  The  oscillations  are  provided  by  a 
current  source  (Ipump) ,  and  may  enter  either  a  shunt 
compliance  provided  by  the  central  airways  (Caw)  or  may 
proceed  further  in  the  axial  direction  toward  the 
alveoli  (whose  compliance  is  Crs) .  Within  this  axial 
path,  inspiratory  flow  is  separated  from  expiratory 
flow  by  idealized  diodes.  in  each  inspiratory  path, 
the  inspiratory  resistance  is  identical  (Rinsp) .  The 
three  models  differ  only  in  the  expiratory  limb  of  the 
axial  path  (enclosed  within  dashed  lines) .  Model  (a) 
includes  a  fixed  expiratory  resistance  (Rexp) ,  which  is 
5  times  larger  than  Rinsp.  Model  (b)  includes  a 
variable  expiratory  resistance  [Rexp (Palv) ] ,  which 
decreases  as  lung  volume  is  increased;  its  value  is  5 
times  that  of  Rinsp  when  Palv  =  +30  cm  H20  total  lung 
capacity  [TLC]  and  10  times  Rinsp  when  Palv  =  0  cm 
H20  (FRC) .  Model  (c)  includes  an  expiratory 
resistance  equal  to  Rinsp,  but  also  contains  a 
transistor  which  models  flow  limitation  so  that  the 
maximum  possible  flow  through  the  expiratory  limb 
increases  linearly  as  lung  volume  (reflected  in  Palv) 
increases.  Values  of  Vcrit  and  the  flow  limiting 
resistor  were  chosen  so  that  Vmax  at  FRC  is  1.26 
liter/sec  and  Vmax  at  TLC  is  2.51  liter/sec. 

Predictions  of  these  models  are  shown  in  Figure  4. 


Figure  3 


Elevation  of  lung  volume  above  resting  FRC  during  HFV 
as  a  function  of  oscillation  frequency,  with  fixed 
tidal  volume.  With  low  mean  airway  pressure  (~  3  cm 
H20)  (triangles  and  solid  lines) ,  there  is  a  steady 
increase  in  dynamic  hyperinflation  with  frequency. 

With  high  mean  airway  pressure  (-  10  cm  H20;  squares 
and  broken  lines) ,  the  pattern  differs:  Although  the 
high  mean  airway  pressure  has  itself  increased  lung 
volume,  no  further  increase  occurs  with  frequency  until 
a  critical  frequency  (2-5  Hz)  is  reached.  As 
oscillatory  frequency  is  increased  further,  dynamic 
hyperinflation  occurs,  in  a  manner  parallel  to  that 
seen  with  low  mean  airway  opening  pressure.  Note  that 
the  critical  frequency  for  dynamic  hyperinflation  at 
high  mean  airway  opening  pressure  is  not  the  same  as 
the  frequency  at  which  Vqq2  plateaus  (see  Figure  1) . 
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Figure  4 

Predicted  dynamic  hyperinflation  (manifest  as  increased 
Palv  above  0  cm  H20  [resting  FRC])  versus  frequency 
at  fixed  tidal  volume  (100  ml)  for  each  of  the  three 
models  shown  in  Figure  3.  Both  models  (a)  and  (b) 
display  progressive  dynamic  hyperinflation  with 
increasing  oscillatory  frequency  at  either  low  mean 
airway  pressure  (0  cm  H20;  squares  and  dashed  lines) 
or  high  mean  airway  pressure  (10  cm  H20;  triangles 
and  solid  lines) .  Only  model  (c) ,  which  incorporates 
expiratory  flow  limitation,  displays  the  plateau  of 
lung  volume  over  the  low  frequency  range  with  high  mean 
airway  pressure,  and  the  similar  dynamic  hyperinflation 
above  a  critical  frequency,  that  were  seen  in  our 
patients  (see  Figure  3) . 


Figure  5 


Branching  network  used  in  the  experiments.  Shown  is 
the  three-generation  network.  The  four-generation 
network  contains  one  additional  generation  with  16 
pistons . 


Figure  6 


Typical  concentration  traces  obtained  by  the  ensemble 
averaging  technique.  Subscripts  refer  to  measurement 
locations  indicated  in  Figure  1.  (a)  Concentration 

traces  for  3  =  1.52,  Pe  =  31.2.  (b)  Concentration 

traces  for  8  =4.81,  Pe=2.35. 
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Figure  7 


Comparison  of  all  experimental  data  to  the  regression 
formula  D  =  0.078  Pe1'66  3-1-45. 


40 


Figure  8 

A  comparison  of  our  results  for  series  A  and  B  to  the 
theoretical  prediction  for  a  straight  tube  with 
sinusoidal  flow  (solid  line).  82%  of  .the  data  fall 
within  the  shaded  zone. 
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Figure  9 


Data  correlations  found  in  experiments  conducted  in 
ranching  networks.  Steady  flow  in  gases  -  solid  line 
(Scherer  et  al,  1975):  oscillatory  flow  in  liquids  - 
dashed  lines  that  bracket  the  range  4  <  S  <  10  (Tarbell 
et  al,  1982)  :  and  oscillatory  flow  in  gases  -  dash-dot: 
line  that  bracket  the  range  4  <  8  <  10  (both  series  A 

and  B  for  this  study) .  Correlations  are  plotted  only 
for  the  range  of  parameters  common  for  all 
experiments.  See  text  for  a  description  of  the 
comparison. 
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